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During the twelve-month period March 1, 1987 through 
February 29, 1988 the focus of activity in the MSMSR program was 
muon spin rotation studies of superconducting materials, in 
particular the new high-Tc materials and the heavy-fermion 
materials CeCu 2 . l Si2 . 

A set of MuSR experiments was done at the Alternating 
Gradient Synchrotron of Brookhaven National Laboratory in May 
1987. Three materials were studied: (1) the high-Tc 
superconductor YBa2Cu307 , (2) an oxygen-depleted variant on the 
"parent" compound to Lai . 8 5 Sro. l 5 Cu04 , La2Cu04-y, and (3) the 
heavy-fermion superconductor CeCu 2 .iSi 2 . 

The data obtained on YBa2Cu307 were combined with data 
obtained on Lai . 8 5 Sro. 1 5 Cu04 in February 1987 in a paper, 
"Systematic Variation of Magnetic-field Penetration Depth in 
High-Tc Superconductors Studied by Muon Spin Relaxation" , which 
has been submitted to Physics Review B for publication. It is 
included as Appendix A to this report. 

The data obtained on La2Cu04-y are presented in a paper 
which was published in Physical Review Letters on August 31, 
1987. This is included as Appendix B to this report. A further 
comparison of these data with neutron scattering results was 
reported by our colleague Y. J. Uemura at the Interlachen HTSC- 
M^S Conference in February 1988. This paper, which will appear 
in the conference proceedings, is also included in this report as 
Appendix C. 

The director spent the period August 1-6, 1987 at the TRIUMF 
laboratory in Vancouver, B.C., working on a MuSR study of 
YBa2 Cu3 0? - y , which was found to be antiferromagnetic below about 
220K, similar to La2Cu04-y. These results were published in the 
March 14, 1988 issue of Physical Review Letters and are included 
in this report as Appendix D. 

The May 1987 study of CeCu 2 .iSi 2 indicated the simultaneous 
presence of superconductivity and static magnetic ordering in 
this material below about 0.8K. Further MuSR studies with a 
helium-dilution refrigerator were clearly needed, so the 
director, along with two colleagues from the College of William 
and Mary (W. J. Kossler and X. H. Yu) and Y. J. Uemura of BNL 
went to the Swiss Institute for Nuclear Research in November 1987 
to study CeCu 2 .iSi 2 in the MuSR low-temperature facility there. 



The phenomenon observed at BNL was confirmed and the magnetic 
ordering appears to be a spin glass. These results are detailed 
in a paper which has been submitted to Physical Review Letters 
and is included here as Appendix E. A shortened and revised 
version of this paper was presented at the Interlaken HTSC-M 2 S 
meeting by Y. J. Uemura. It is included here as Appendix F. The 
travel to Switzerland was supported by the National Science 
Foundation. 

Apart from superconductors, we also continued to analyze 
data taken during earlier experimental runs on metal hydrides. A 
paper on "Muon Motion in Titanium Hydride" is being submitted to 
Physical Review B and is included here as Appendix G. 

The director spent the periods July 11-August 1 and August 
12-September 3, 1987 at the Los Alamos Meson Physics Facility 
(LAMPF) working on two nuclear physics experiments, a search for 
eta-mesic nuclei, and a study of nuclear pion charge -exchange 
reactions at rest. His travel was supported by a National 
Science Foundation grant administered by the College of William 
and Mary. VSU graduate student Norman Fuqua (since graduated) 
participated in the July run. 

A paper based upon a December 1986 search for eta-mesic 
nuclei done at the Brookhaven AGS has been accepted for 
publication in Physical Review Letters. It is included in this 
report as Appendix H. 

A long set of MuSR experiments at BNL began in February 1988 
and continued into May 1988 (which explains the lateness of this 
annual report). A variety of high-Tc materials were studied, 
including YBa 2 Cu3 - x Cox 07 . These data will be described in the 
next six-month report. 

Lucian R. Goode, Jr. , who did his thesis research on MuSR 
studies of iron alloyed with gold, dysprosium and tantalum, 
received the M. S. degree in physics in May 1987. Norman Fuqua 
wrote his thesis on inelastic scattering of polarized protons 
from 12 C (done at Los Alamos) and received the M. S. degree in 
July 1987. Li-Tai Song has passed the orals on his thesis, "The 
Surface Muon Beam at Brookhaven National Laboratory" and should 
receive the M. S. degree in July 1988. Nana Adu has passed the 
orals on his thesis, "A Muon Spin Rotation Study of Zirconium 
Hydride", and should also receive the M. S. degree in July 1988. 
A senior physics major, Michael R. Davis, who worked extensively 
on this research program, graduated magna cum laude from VSU in 
July 1987. 

The director gave four talks on the studies of 
superconducting materials done at BNL, TRIUMF, and SIN during 
this reporting period. They were presented at the Virginia 
Academy of Science meeting in Norfolk in May 1987, at a MuSR 
colloquium held at the Los Alamos Meson Physics Facility in 
August 1987, at a NASA workshop held at the NASA Langley Research 



Center in January 1988, and at a meeting of the VSU Society of 
Physics Students, also held in January 1988. A presentation was 
planned for the Virginia Academy of Science meeting to be held in 
Charlottesville in May 1988. 

The director continues to serve on the board of trustees of 
the Southeastern Universities Research Association, which is 
building the Continuous Electron Beam Accelerator Facility in 
Newport News, VA, and which operates the SURANet computer 
network. In December 1987 he was appointed to the board of 
directors of the Virginia Research Network, and in May 1987 he 
was appointed to the Governor’s Ad Hoc Committee on the 
Superconducting Super Collider. 

A temperature controller and a CAMAC crate controller, both 
ordered during the SSPRI grant period, arrived during the summer 
of 1987 and have been installed in the MuSR rig at BNL. An 8-MHz 
80287 math chip for the IBM PC/XT has been installed and is 
operating. 

Dr. Robert I. Grynszpan, the co-investigator of this 
research program, was in residence at VSU between October 28 and 
November 14, 1987. He worked with Dr. Min Namkung of NASA 
Langley on magnetoacoustic measurements of fatigue in metals 
during his visit. These results will be included in the next 
report. A student of Dr. Grynszpan from the University of Paris, 
Patrick Langlois, is expected to visit VSU as an exchange student 
during the summer of 1988. He will work with Dr. Namkung at NASA 
Langley on further magnetoacoustic measurements of fatigue in 
metals . 

The report period of March 1, 1987 to February 29, 1988 saw 
the development of muon spin rotation into perhaps the best 
technique for measuring the microscopic properties of high-Tc and 
heavy-fermion superconductors. The coming year is full of 
promise for continued studies in these areas, as well as the more 
traditional areas which utilize the muon probe, such as metal 
hydrides and magnetism. 

We appreciate the support of NASA, and the collaboration of 
scientists from William and' Mary, Brookhaven National Laboratory, 
George Mason University, and NASA Langley, all of which have been 
vital to the success of this program. 


Carey E. Stronach 

Director 

June 3, 1988 
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ABSTRACT 

The muon spin relaxation rate o has been measured in the high-Tc superconductors 
Y BaiCuzOx for z = 6.66, 6.95, 7.0, and Lai.gsSro.isC'uCU in transverse external mag- 
netic fields 1 ~ 4 kG. We find a simple relation which connects the transition temperature 
Tc , the magnetic-field penetration depth A l, the carrier concentration n 3 , and the effective 
mass m* as Tc a c a 1/A l 2 cc n 3 /m*. The linear dependence Tc <x n 3 /m* suggests a 
high energy scale for the coupling between superconducting carriers. 



The discovery 1,2 of the layered Oxide high-Tc superconductor systems La 2 - y Sr y CuO 4 
and Y Ba-iCuzOj has triggered extensive experimental activity 3 and renewed theoreti- 
cal interest 4,5 in the search for a novel mechanism for superconductivity. Muon . Spin 
Relaxation (/iSR) is a direct method 6 for measuring magnetic-field penetration depth 
in superconductors 7 . /xSR has been applied to the high-Tc oxide superconductors 8 " 11 
and related antiferromagnets 12 ’ 13 . Because of recent technological development in sam- 
ple preparation, it has now become possible 14 to study single-phase specimens with the 

oxygen concentration controlled to within ± 0.02 per formula unit (/ f.u.). In this paper, 

«» 

we present /z SR measurements on YBazCuzOx superconductors with averaged oxygen 
concentrations x = 7.0, 6.95, and 6.66 /f.u.. We combine these results with the earlier 
work 9 on Lai.ssSr 0 . i^CuO^ and focus on the systematic dependence of the observed 
muon relaxation rate a and the derived penetration depth A^. The results indicate that 
the superconducting transition temperature Tc is approximately proportional to the su- 
perconducting carrier concentration n 3 divided by the effective mass m * . We discuss the 
implication of this relation on energy scales of the coupling between the carriers. 

The sintered pellet specimen of Y BazCuzO-j was prepared using a method described 
in ref. 15. The powder specimens of YBclzCuzOg.qs and YBa 2 Cu 3 O 6 .se were prepared by 
-using another method described in ref. 14. which reports Meisner effect and susceptibility 
measurements on a series of Y BazCuzOx specimens ranging from x = 6.0 to 7.0. The spec- 
imens with x > 6.5 show superconductivity. The /iSR experiments on Y BazCuzOj were 
performed at the AGS muon channel of Brookhaven National Laboratory with a trans- 
verse external magnetic field H ex t of 1 kG applied perpendicular to the initial direction' of 
muon spin polarization. The measurements on YBazCuzOs.^s and Y Ba 2 Cu 3 O 6 .se were 
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carried out at the M15 muon channel of TRIUMF (Vancouver) with H ex t = 4 kG. In both 
cases, the data were taken by cooling the specimen in external field from T > Tc to lower 
temperatures. 

In the transverse-field fj. + SR experiments, one observes the decay time histogram of 
positive muons stopped in the specimen 

N(t) <x expl-t/r^l + AG x (t)cos(u/ fJL t)\, (l) 

where is the muon lifetime 2.2 nsec, A is the initial precession asymmetry, u/ M is the muon 
precession frequency, and the relaxation function G x (t ) represents the time evolution of the 
muon spin polarization. At all measuring temperatures (5 . OK < T< 300#), the observed 
precession amplitude A indicates that within experimental error all the muons stopped in 
the specimen contribute to the precession signal. The frequency was approximately 
equal to ("V = 2r x 1.355 x 10 4 /Oe is the gyromagnetic ratio of n + ) above Tc- 

ojfj. decreased slightly with decreasing temperature below Tc, due to the partial exclusion 
of the external field H ex t in the type II superconductors at ff e it > #ci* For simplicity, 
here we assume a Gaussian shape for G x (t): 

2 2 

G x {t) = exp(-^~), (2) 

where a is the muon spin relaxation rate. 

Figure 1 shows the temperature dependence of a obtained for the present Y Ba 2 Cu^O x 
compounds together with the earlier results 9 on Lai.8s5r 0 .isC , uO4. The very small values 
of a observed in all the specimens above Tc can be accounted for by nuclear dipolar 
broadening. Combining this feature with the full amplitude for A and the reduction 



of Up below Tc, one can conclude that there is no static magnetic ordering in these 
superconducting specimens either above or below Tc- This aspect was comfirmed in the 
zero-field /x SR measurements on Y BazCuzO? . Below Tc, the value of a increases rapidly 
with decreasing temperature. This is due to the inhomogeneity of the static local field at 
the muon site in the type-II superconducting state where H ext penetrates as a lattice of 
flux vortices. We notice here that the four different specimens in Fig. 1 have reasonably 
similar shapes for the curvature of the temperature dependences cr[T). This implies that 
Tc is approximately proportional to o(T — > 0), as demonstrated in Fig. 2 for the four 
different specimens. 

Pincus et al. 16 used the London equation to calculate the distribution of magnetic 
fields in the vortex state, and obtained the second moment 


yj {AH 2 ) = 


00 

A|\/l67r 3 


( 3 ) 


with the flux quanta 0 q, for the square lattice of the vortex when the second moment 
becomes independent of the external field H ex t , i.e., when \l is comparable to or greater 
than the distance between adjacent vortices. The present condition with H ex t = 1 ~ 4kG 
satisfies this criterion. For a triangular lattice, one needs to multiply 10 a factor of 0.93 
to the right hand side of eq. (3). Then, one can deduce the value of the penetration 
depth directly from the observed relaxation rate a which corresponds to ~\p\J (AH 2 ). 
Figure 3 shows the temperature dependence of thus obtained for the triangular vortex 
lattice. The values for A i(T — ► 0) are- listed in Table 1. In the field-cooled measurements, 
the density of the magnetic flux is kept almost constant above and below Tc- If one 
changes the external field in the superconducting state, in contrast, the flux vortices have 
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to move within the sample to change the spatial flux density, and thus the experimental 
results become sensitive to the flux pinning 17 . Therefore, it is important to measure the 
penetration depth in the field-cooled condition, as in the present experiment. 

In actual systems, we noticed that the functional form of G x (t ) is somewhat in between 
Gaussian and exponential. This is due to the complicated distribution of magnetic fields 
for the vortex lattice as well as to the effect of anisotropy on the penetration depth Ax,. 
Correcting for the former effect would reduce the resulting values of Ax, by about 30 %. For 
the case of maximum anisotropy where the penetration depth A in the soft direction ( H ex t 
applied parallel to the CuO plane) is infinite, Celio et al. 18 find that the value of A for 
H ex t J_ {CuO plane} to be about half the powder averaged value. These corrections make 
it difficult to deduce the absolute values of A*, accurately. However, we stress here that 
the relation a oc l/A^ 2 holds for any of the above calculations. Therefore, the systematic 
and temperature variations of Ax, can be discussed based on Figs. 1 and 3. 

The London penetration depth Ax, is given as a function of effective mass m* and the 
carrier density n s as 


A r = 


m*c 2 
47 t n,e 2 ’ 


( 4 ) 


Combining this with the relation cr = 7 M -y/(Aif 2 ), eq.(3), and the approximate experi- 
mental result <r a Tc, we obtain a simple- relation 


a cc 


^ rp 

T y °c — a T c . 
A l rn* 


( 5 ) 


Thus, the transition temperatures Tc of the four different samples are simply propor- 
tional to the carrier concentration n 3 divided by the effective mass m* regardless of the 
crystallographic differences of the samples. 
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The relation. Tc oc n (n: normal-state carrier density) has been found by Hall con- 


stant and related measurements for La 2 - y Sr y CuO 4 between y‘ = 0 ~ 0.15 (ref. 19) and 
Y BaiCuzOx between x = 6.5 ~ 7.0 (ref. 20). This linear relation can also be obtained 
in a calculation of the number of holes for a formula unit, assuming charge neutrality for 
La. 2 Cv.O 4 and YBa 2 Cu z O QA and adding 1 hole for the substitution of Sr to La and 2 
holes for additional oxygen per formula unit. These results suggest that the carrier con- 
centration n a , rather than the effective mass m*, plays a major role in changing Tc in eq. 
5. Indeed, the Pauli susceptibility at T > Tc or the Sommerfeld constant 7 of the low 
temperature specific heat, which are proportional to m’n 1 / 3 , do not depend much on the 


differences in material 21 , supporting the above view point. The present work has shown 
that the linear relation holds for the superconducting carrier concentration n 3 . 

We now discuss the implications of the relation Tc oc n a /m*. In the BCS theory 22 
with the phonon-mediated coupling of electrons, Tc is given in the so-called weak coupling 
limit as 


fc bTc 


2 hu>D ezp( 




( 6 ) 


where kg is the Boltzmann constant, u>g is the Debye frequency, V represents the effective 
attractive interaction, and.J 9(e/) is the density of states at the Fermi energy e/. To obtain 
this equation, one assumes e/ kug and solves the secular equation by integrating the 
energy range of the coupling interaction 0 — * hup, which results in the pre-exponential 
factor huig- It is difficult to reconcile the relation Tc oc n s /m * with eq. ( 6 ), because 
| D(e/) does not depend on n in the 2-dimensional non-interacting electron gas. Recent 
single-crystal measurements 23 on H c 2 and on the transport properties suggest a highly 
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2 -dimensional character for the electron system. Moreover, the magnitude of the electron- 
phonon interaction V, inferred from the temperature dependence of linear resistivity in 
the normal state 24 , is too small to explain the high transition temperature in the standard 
phonon-mediated mechanism. 


In contrast, when the energy scale of the attractive interaction which couples the 
carriers is larger than that of ep, the energy integration in the secular equation runs over 
the range 0 to e/. This would, put, roughly speaking, £/ in the pre-exponential factor 
of eq. (6) instead of hup. In a two-dimensional. non-interacting electron gas, the Fermi 
energy €/ is proportional to the quantity n/m* . Then one could expect the simple relation 
Tc oc n s /m*. This argument works without essential change also for three dimensional 
systems where ep cc n 2 / 3 /m*. Thus, the relation Tc oc n 9 /m * suggests a high energy 
scale of the interaction which mediates the coupling between superconducting carriers in 
high- Tc superconductors. Such a high energy scale may be found in models based on the 
large transfer integral of a carrier between the oxygen and neighboring copper atoms 5 . 
The linear relation Tc oc n 3 is also expected in a resonating- valence-bond picture 25 . 

We study the temperature dependence of \i with examples of YBa^Cu^Os ^ 5 and 
Y Bo.2Cu2O7.Q~ The sharp changes of o{T) near Tc observed for these samples suggest a 
good homogeneity in the oxygen concentration. As shown by the solid lines in Fig. 3 , the 
experimental data agree well with the empirical formula 22 


A l(T) 


Al(T = 0) 

^l-P/Tc)*' 


(?) 


This result is consistent with earlier y SR. works 8,10 , but disagrees with a recent bulk 
measurement 26 . Equation ( 7 ) is calculated for much smaller than the coherence length 
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£ by assuming an isotropic energy gap at the Fermi surface. For A l » the BCS theory 22 
predicts that Ax,(T) increases more rapidly than eq. (7) with increasing temperature at 
T < OJTc • The anomalous zeros of the energy gap at some point or line of the Fermi 
surface would change the theoretical curves for A l(T) to increase faster with increasing 
temperature in the low temperature region. Therefore, the present results suggest that 
the energy gap is predominantly finite. 

In summary, based on the muon spin relaxation experiments, we have shown that the 
approximate proportionality Tc ct n s /m * holds universally for different high-Tc oxide su- 
perconductor systems. This feature suggests the high energy scale of the coupling between 
superconducting carriers. 

This work is supported by the Division of Materials Sciences, US Department of 
Energy under contract 76-AC02-CH00016, the US National Science Foundation under 
DMR 8503223, NASA under NAG- 1-416, and by the NSERC of Canada. 
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FIGURE CAPTIONS 


Fig. 1. 

Temperature dependence of the relaxation rate a of the muon spin polarization, as 
defined in eq. (2), observed in four different specimens of high - Tc superconductors. Data 
on LaiAsSrQ'isCuOt and Y Ba. 2 Cv.zO 7 were obtained in a transverse external magnetic 
field of 1 A:C?, while the measurements on the other two specimens were performed with a 
field of 4 kG. 


Fig. 2 

Superconducting transition temperature Tc, as determined by the /xSR measurement, 
plotted versus the values of the muon relaxation rate a at T — + 0 for the four different 
specimens of high - Tc superconductors. 


Fig. 3 

Temperature dependence of the magnetic field penetration depth A 1 derived from the 
muon spin relaxation rate shown in Fig. 1. was calculated using a simple approximation 
for the triangular vortex lattice as described in the text. The solid lines represent fits of 
the data to eq. 7, with A i{T = 0) = 1656 A and Tc — 89.9 K for Y Ba. 2 Cv. 2 OQ as, and 
A £,(T = 0) = 1472 A and Tc — 91.1 K for Y Ba. 2 Cv. 2 O 7 A. 
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ABSTRACT 


Muon spin rotation and neutron scattering studies on powder and single- crystal specimens 
of La^CuO^-y are compared. The apparent difference between the muon and neutron 
results for the ordered moment in the antiferromagnetic state is interpreted as the signature 
of increasingly short-ranged spatial spin correlations with increasing oxygen content. 

(submitted to the Interlaken ETSC-M^S Conference, February, 1988) 
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It has been known 1 that the magnetic properties of the antiferromagnetic compound 
La. 2 Cv.O 4 -y depend sensitively on small differences in the oxygen content y = 0 ~ 0.03. 
The input of 0 2 “, equivalent to the substitution of Sr 2+ for ia 3+ , removes electrons and 
creates holes in the system, and suppresses the 3-dimensional antiferromagnetic ordering. 
To study this phenomenon with direct microscopic magnetic probes, experiments were 
performed using the muon spin rotation (pSR) (ref. 2) and neutron scattering 3 techniques. 
In this paper we present new jiSR data obtained at TRIUMF on a single-crystal and several 
powder specimens, compare the muon and neutron results, and discuss their implications. 

Figure 1 shows the muon spin precession frequency i/ M observed in zero field for four 
different powder specimens of La.2Cv.O4-y. i/ M is proportional to the static internal mag- 
netic field at the muon site from nearby ordered magnetic moments. Therefore, Fig. 1 
represents the temperature and y variation of the ordered Cu moment 5. It is remarkable 
that i / m (T — *■ 0) oc S(T — ► 0) changes only within about 15 % for the four specimens 
despite the large difference in the ordering temperatures Tn (from 300 K to 15 K). The 
ordered moment has also been measured by neutron scattering by Yamada et al. (ref. 3) 
on a few single crystal specimens of La. 2 Cv.O 4 -y. Figure 2 shows the temperature vari- 
ation of the sublattice magnetization M 9 derived from the intensity of the 3-dimensional 
magnetic Bragg reflection. In contrast to the /iSR results, M t [T —* 0) decreases by more 
than a factor of 3 for the decreasing Neel temperatures of the specimens. In this respect, 
H SR and neutron measurements look inconsistent. 

To study further details, we performed fJ.SH measurements on a single-crystal specimen 
of La 2 Cv.O 4 -y (NTT No. 3 specimen in ref. 4, grown simultaneously in the same batch 
as sample No. 3 in Fig. 2 and in ref. 2). The neutron magnetic Bragg intensity was 



measured on this crystal by Endoh zt al . 4 , and the temperature variation was very close 
to the results shown for sample No. 3" in Fig. 2. The ordered moment estimated from 
the neutron study was about Ms{T — ► 0) ~ 0.15/ifl. Figure 3-(a) shows the muon spin 
precession frequency measured on this specimen. i/ m (T — *• 0) is approximately equal to 
the values observed for the other powder specimens (Fig. 1), thus indicating that S(T — ► 0) 
is not much different from the value 0.5 ~ 0.6 he determined by the neutron experiment 
(Fig. 2) on the specimens with TV ~ 300iir. 

We also performed /iSR measurements in a weak transverse magnetic field of about 
100 G, and determined the volume fraction of the magnetically ordered part of the specimen 
as shown in Fig. 3-(b). The magnetic ordering takes place gradually between 200if and 
10QK, below which almost the entire volume of the specimen becomes antiferromagnetic. 
Lack of a sharp magnetic ordering is assumed to be due to the inhomogeneous spread of 
the oxygen concentration in the large single crystal specimen. Similarly, it is confirmed 
that the predominant volume fraction of the powder specimens orders magnetically below 
TV. This rules out the explanation that the difference between muon and neutron results 
may be due to a non-order ed volume fraction of the specimens. 

The difference between /iSR and neutron results can be understood in the following 
way. The /x + is a point-like probe in real space, and the local field at a muon site is 
due mostly to its neighbouring Cu moments. Therefore, once the antiferromagnetic spin 
configuration becomes static ( t > 1/zsec), even if the spatial correlation may be short 
ranged, the /i + sees the specimen to be almost perfectly antiferromagnetic. In contrast, 
the elastic magnetic scattering intensity of neutrons forms a peak in reciprocal space at the 
Bragg point only when the spatial 3-dimensional ordering is long ranged. Therefore, the 


present results suggest that the magnetic ordering becomes more and more short-ranged 
with increasing oxygen content (i.e., with increasing number of holes), resulting in the 
decreasing Neel temperatures. 

The closed-shell electron configuration of 0 2 ~ mediates the super-exchange antiferro- 
magnetic coupling between adjacent Cxi moments. When the introduction of a hole makes 
the oxygen to be O x ~ (ref 5), the unpaired spin at the oxygen would mediate the effec- 
tively ferromagnetic coupling between neighbouring Cxi moments. Therefore, the input of 
oxygen would create frustrated bonds 6 in the Cxi — 0 plane, and thus help to destroy the 
long-ranged spin correlation. This picture provides a possible explanation for the present 
results. There are two possiblities for the short-ranged correlation: a) the randomness 
introduced within the Cxi — O plane; and b) the cut-off of the spin correlation between 
different Cxi — O planes. It is not possible to determine from the present work whether 
either or both of them take place. The neutron study of ref. 4 gives support to b) for the 
particular case of the specimen shown in Fig. 3. 

By combining neutron and muon experiments, one can estimate the microscopic static 
moment to be S = 0.5 ~ 0.6^s /Cu. This value is significantly smaller than the integer 
moment for spin 1/2, but it can be explained by the quantum spin reduction for the spin 
1/2 2-d Heisenberg system obtained from the statistical-mechanics calculations 7 . A spin- 
density-wave state is another possibility to create the non-integer moment. From the 
present work alone, it is not posssible to tell which is the case. 

We are grateful for useful discussions with R.J. Birgeneau, V.J. Emery, G. Shirane 
and S.K. Sinha. This work is supported by the USDOE (DE-AC02-76CH00016), the US 
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NSP (DMR 8503223), NASA (NAG- 1-416), the Canadian NSERC, and by the Japanese 
Ministry of Education, Culture and Science (Grant-in- Aid for Scientific Research). 
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Figure Captions 


Fig. 1 Muon spin precession frequency observed in zero external field on various sintered- 
powder specimens of La 2 CuO^- y . 

Fig. 2 Sublattice magnetization M a derived from the magnetic Bragg-peak intensity of 
neutrons on single-crystal specimens of La^CuO^y by Yamada et al. (after ref. 3). 


Fig. 
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(a) Muon spin precession frequency i/ M measured in zero external field on the NTT 
No. 3 single-crystal of La^CuO^y. (b) Volume fraction of the paramagnetic region 
of the specimen measured by muon spin rotation experiments in a weak transverse 
external magnetic field. It is shown that almost the entire volume of the sample orders 
magnetically below about T ~ 100 JC. 
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ABSTRACT 

Zero- and longitudinal-field muon spin relaxation measurements on a heavy-fermion system 
CeCu<z.iSi 2 have revealed an onset of static magnetic ordering below Tm ~ 0.8FC, which 
coexists with superconductivity below Tc = 0.7 K. The line shapes of the observed muon 
spin depolarization functions suggest an ordering in either spin glass or incommensurate 
spin-density-wave state, with a small averaged static moment of the order of 0.1/zb per 
formula unit at T -» 0. ' (PACS Nos. 75.20.H, 76.75, 74.70.H) 
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C zCuiSi-i belongs to a group of Cerium or Uranium intermetallic compo un ds, called 
heavy-fermion systems 1 , which are characterized by the extremely large T-linear term 
C = 7T with 7 ~ 1 Joule/ mole deg 2 of the electronic specific heat C at low temperatures 
inferring very large effective mass m*. In 1979, Steglich et al 2 found a superconducting 
transition in C eCu^Sii at T ~ 0.5if , as the first example of the superconducting ground 
state of a heavy-fermion system. This discovery triggered extensive study of supercon- 
ductivity and its possible relation with magnetism in the highly correlated heavy electron 
systems 1,s . Although there are some susceptibility measurements 4,5 which suggest mag-' 
netic ordering of non-superconducting systems C ei^yLayCu^Si^ with y > 0.2 (ref. 5) 
and CeCu x Si2 with x = 1.9 (ref. 4), so far it has been common to assume 1 a purely 
superconducting ground state without magnetic ordering for the superconducting com- 
pounds CeCu x S i2 with 2.0 < x < 2.2 (ref. 6). In this paper, we present the first direct 
evidence from zero-field muon spin relaxation measurements that superconductivity and 
static magnetic ordering coexist in CeCu2.iSt2 • Muon Spin Relaxation (//SR) 7 is a very 
powerful tool to detect static magnetic ordering 8,9 . Previous applications of //SR to heavy- 
fermion systems led to discoveries of magnetic ordering with quite small ordered moments 
(0.001 ~ 0.05 iib) in superconducting UPtz (ref. 10) and non-superconducting CeAlz (ref. 
11). With the present results on CeCu2.-i.Si2, we now have three heavy-fermion super- 
conductors, CeCu2.iSt2,UPt3, and URU2S12 (ref. 12), which show co-existing magnetic 
orderings. 

Superconductivity of the stoichiometric CeCu2.oSt2 is known to be somewhat unstable 5 . 
A small amount of additional Cu helps to stabilize the superconductivity of CeCu x Si2 
with x = 2.1 ~ 2.2 whose superconducting transition temperature Tc is around 0.7 K 
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(ref. 6). We therefore prepared a poly-crystalline sample of CeCu 2 .iSi 2 by the method 
described in ref. 6. Figure 1 shows resistivity measured on a piece cut off from the sample 
of CeCu 2 .\Si 2 used in the present /z SR measurement. Reference 6 describes detailed resis- 
tivity and specific-heat measurements on a few poly-crystalline samples of CeCu x Si 2 with 
x = 1.9 ~ 2.2 made with the same method. The present specimen has Tc at around 0.7 K , 
and the temperature dependence of resistivity shows the same curvature as reported in ref. 
6. To further characterize the present sample, we made a neutron scattering measurement 
of the crystal structure, and confirmed that it is a single phase material without any minor 
phase within the experimental accuracy of a few volume percent. 

We started zero-field /zSR measurements on the present sample of CeCu 2 .iSi 2 at the 
AGS muon channel of Brookhaven National Laboratory by using a 3 He- cryostat. A rapid 
increase of the muon spin relaxtion rate was found with decreasing temperature below 
T ~ 0.8i£\ We then continued the measurement at SIN, Switzerland, by using a surface 
muon beam and a dilution refrigerator. The facility at SIN allowed full access to a wide 
temperature region with high-statistics data as we report in this paper. A positive muon 
beam was stopped at the sample of C eCu 2 ,\Si 2 [2cm x 1cm x 0.5cm) mounted on a cold 
finger of the dilution cryostat, and the muon-decay positrons were recorded mainly with a 
counter placed in the forward direction with respect to the beam direction. The counting 
rate I(t) of this counter is given as 

I(t) cc exp(-t/r M )[l - AG z (t)\, (1) 

where is the muon lifetime 2.2 /zsec, A is the initial decay asymmetry (A ~ 0.25- at 
the present condition), and the relaxation function G z [t) represents the time evolution of 
muon spin polarization. 


3 



Figure 2-(a) shows the relaxation function G z (t) thus observed in zero field. The 
relaxation rate increased rapidly below T ~ 0.8 K with decreasing temperature. In general, 
the depolarization of muon spins in zero field can be due either to randomness of the static 
internal local field Hint or to the fluctuating dynamic local fields. One can distinguish 
between these two cases by making measurements in the longitudinal external magnetic 
field H ex t applied parallel to the initial muon spin direction (i.e., the direction of muon 
beam, denoted as z direction hereafter). When the internal field is static and H cz t is 
larger than Hint, H ex t can align the local field Hi oc = Hint + H ex t to be nearly parallel 
to the muon spin polarization, thus keeping G z (t) finite. In contrast, the dynamic spin 

fluctuations are usually much faster than the corresponding Zeeman frequency ui = 7^ H ex t 

\ 

(7^ = 2 tt x 13.554MHz/kG) of fx + , so that there is almost no effect of H ext on G z (t) for 
the dynamic case. We have performed such measurements in H ex t = 250G and IkG at 
T = 0.1 K as shown in Fig. 2~(b). The longitudinal fields suppress the depolarization 
and change G z {t) remarkably. This indicates that the depolarization observed in zero 
field is due mainly to the static random local fields of the order of 100 ~ 200 G. Similar 
measurements with H ex t confirmed that the depolarization at T = 0.8 K is also due to the 
static fields. 

The line shapes of G z (t) in zero field shown in Fig. 2-(a) resemble those observed 
in the dilute-alloy spin glasses CuMn or .AuFe (ref. 8). The lack of coherent oscillation 
indicates that the magnitude of the static internal field | Hint | has a wide distribution. 
For a uniform | Hint |, one should have seen a muon spin precession as reported in refs. 9 
and 11. G z (t) in Fig. 2-(a) at T = 0.05IiC starts with a Gaussian-like shape at t — + 0. This 
is somewhat different from the case in the dilute-alloy spin glasses where G z (t —* 0) decays 
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with an exponentiai-like shape when the temperature is well below the susceptibility-cusp 
temperature T g . In this paper, we do not develop a complicated theory to account for this 
line shape, but rather adopt a phenomenological relaxation function 

Gz{t) = ^exp(-i<7 2 t 2 ) + ^exp(-At) (2) 

to fit the observed data. The first term corresponds to the quick initial decay caused by 
the x- and y- components of the static random local fields. The second term represents 
the tail arising from the component of U, nt parallel to the initial muon spin direction 
(i.e., the 2 -direction). Reasonably good fits to the data at all different temperatures 
0.05if < T < l.Oif were obtained when we assumed A\/A = 2/3, A^/A = 1/3, with 
A = 0.25. The present counter configuration without the backward counter, however, 
made it difficult to determine the second term of eq. (2) accurately. Therefore, the decay 
rate A of the tail, ranging around A ~ l.Opsec -1 , may be subject to a large systematic 
error. In contrast, the first term of eq. (2) can be determined with reasonable precision. 
From the fit to eq. (2) represented by the solid lines in Fig. 2- (a), we thus obtained the 
muon spin relaxation rate o as shown in Fig. 3. 

The relaxation rate o in Fig. 3 increases rapidly with decreasing temperature below 

j 

T <— 0.8K. This indicates the onset of a random but static magnetic order at the magnetic 
ordering temeprature Tm ~ 0.8K. Because of the limited accuracy of the temperature 
measurements and control with the cold-finger cryostat, as represented by the horizontal 
error-bars in Fig. 3, we can not identify whether the magnetic and superconducting or- 
derings occur simultaneously or whether Tm and Tc are somewhat different. It is evident 
from Figs. 1-3, however, that the superconductivity and magnetic ordering coexist below 
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T ~ 0.7 K in CeCu 2 .iSt 2 . The present results indicate that more than 90 vol um e percent 
of the specimen undergoes the magnetic ordering. The relaxation rate a in Fig. 3 repre- 
sents a measure of the width of the static internal fields as a ~ 7ft \/< (A Hint) 2 >. The 
observed value of o ~ \0jMsec ~ 1 at T — ► 0 corresponds to the static random local fields 
of the order of o/ / 7 #t ~ 120G. In most magnetic materials, local fields at muon sites axe 
due mainly to the dipolax field from the surrounding moments. Since ) a + is a point-like 
probe in real space, and since information on the muon stopping site is lacking, it is not 
easy to accurately estimate the spatial spin structure of the magnetically ordered state of 
CeCu 2 .\Si 2 . It is, however, possible to note the following. 

Spin glass ordering is one of the most likely spin structures which produce a laxge 
distribution of Hint as observed in the present experiment (see ref. 8). This picture is 
consistent with the observation of a strong effect of external fields on the susceptibility 
measurements of CeCui,gSi 2 (ref. 4). If one assumes that the majority of Ce (or Cu ) 
atoms participate in the magnetic ordering, the observed width — 120G of Hint corresponds 
to the dipolax field from the ordered moment of the order of 0.1 /zb . If instead the small 
population of Ce 3+ ions forms a spin glass with an ordered moment of 5 /XjS> like dilute- 
alloy spin glasses, then the observed value of a corresponds to an ordering of about a 
few percent of the Ce atoms. The latter picture is advocated in ref. 4 for the case of 
CeCu\.gSi 2 . The Gaussian-like decay of G z (t — ► 0) may favor the former type of spin 
glass, but then one faces the difficult question as to the origin of the frustration of the 
exchange interactions. An incommensurate spin-density-wave state, like that observed in 
CePbz (ref. 13), is another possible spin structure which gives the local field at the muon 
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site a wide distribution. Unfortunately, the present experiment alone can not distinguish 
the above mentioned three possible spin structures. 

A neutron scattering experiment on CeC« 2 .iSt 2 is underway 14 to study spatial spin 
correlation. For such a small averaged moment as ~ 0.1 hb , however, the magnetic scatter- 
ing intensity of neutrons is very small. Neutron measurements become even more difficult 
when the spin glass ordering makes the scattering diffusive in reciprocal space. In this 
respect, the present experiment demonstrates the unique capability of zero-field /J.SR to 
detect magnetic orderings with small averaged moments. In the previous /zSR study on the 
same specimen of C eCu 2 .\St 2 performed in the transverse external magnetic field 15 , the 
depolarization of muon spins observed below T ~ 0.8if was tentatively attributed to the 
inhomogeneous penetration of the transverse external field in the type- II superconducting 
state. Zero-field n SR can detect the magnetic ordering without the complication of the 
field penetration. The present results indicate that a major part of the depolarization 
observed in the transverse field was due to the static magnetic ordering. In order to find 
out whether the superconductivity and magnetic ordering occur at the same tempeature 
or at different temperatures, we are planning to perform additional /zSR measurements on 
specimens with different Cu stoichiometry which may have different Tm (see ref. 4). 

Recently, a heavy-fermion superconductor UPtz ( Tc ~ 0.5 K) was found 10 to order 
magnetically below T ~ 5 K with an extremely small averaged moment of ~ O.OOI^b- 
URu 2 S %2 is another superconducting heavy-fermion system (Tc ~ l.Oif) which orders 
antiferromagnetically below Tu <— 17 K with an ordered moment of 0.03 per Uranium 
atom 12 . TJBeiz is so far the only remaining heavy-fermion superconductor without an 
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identified magnetic ground state, yet there is a report 16 which suggests a possible antifer- 
romagnetic ordering in UBt\z doped with a small amount of Th. Together with the present 
results on CtCu 2 .\S%i , these features indicate that the magnetic ordering with extremely 
small averaged moment may be a common feature of heavy-fermion superconductors. 

In summary, zero-field y SR measurements on CeCu 2 .\Si 2 have shown clear evidence 
of static magnetic ordering below Tm ~ 0.8if with a very small averaged moment of 
~ 0.1 he in either a spin glass or an incommensurate spin-density-wavestate. This ordering 
coexists with superconductivity below Tc = 0.7JFC. Further experimental and theoretical 
studies are clearly needed for the full understanding of the role of such magnetic orderings 
on the superconductivity of heavy-fermion systems. 

We would like to acknowledge useful discusions with P.B. Allen, K. Kakurai, M. Steiner 
and E. Recknagel. We thank W. Schonig for help on the measurements at Brookhaven. 
This work is supported by the Division of Materials Sciences, US Department of Energy 
under contract 76-AC02-CH00016, the National Science Foundation under DMR 8503223 
and INT 8413978, NASA under NAG-1-416, and by the Japanese Ministry of Education, 
Science and Culture under the Grant-In-Aid for Scientific Research. 
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FIGURE CAPTIONS 


Fig. 1. Resistivity measured on a piece of CeCu^.xSi^ cut off from the present specimen used 
in the fiSH measurements. Superconducting transition occurs at Tc = 0.7FC . 

Fig. 2. (a) Muon spin relaxation function G z [t) in zero field observed in CeCui.iSii. Solid 
lines represent fits to eq. (2). (b) Muon spin relaxation function in CeCu 2 .xSt 2 

observed at T = 0.1 K in longitudinal external magnetic fields LF of 0G, 250G, and 
IkG. Solid lines are guides to the eye. 

Fig. 3. Muon spin depolarization rate a, as defined in eq. (2), derived from the relaxation 
functions observed in CeCu 2 .iSi 2 in zero field. The onset of magnetic ordering is seen 
around Tm ~ 0.8JT. 
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ABSTRACT 

Zero- and longitudinal-field muon spin relaxation measurements on a poly-crystal sample 
of a heavy-fermion superconductor CeCu 2 .iSi 2 ( Tq = 0.7 K) have revealed an onset of 
static magnetic ordering below T ~ 0.8K. The line shapes of the observed spectra in zero 
field indicate a wide distribution of static random local fields at muon sites, suggesting 
that the ordering is either spin glass or incommensurate spin-density-wave state. The 
observed width of the random local field at T = 0-05J5C corresponds to a small averaged 
static moment of the order of Q.lps per formula unit. 




C eC u^S t *2 is the first heavy-fermion system which was found to become superconduct- 
ing at Tc ~ Q.5K (ref. 1). Although there were signatures suggesting possible magnetic 
orderings in non-superconducting CeCui.gS;^ (ref. 2) and Cei- y La y Cu 2 St 2 with y > 0.2 
(ref. 3), superconducting specimens CtCu x S \2 with x — 2.0 ~ 2.2 have so far been be- 
lieved to have purely superconducting ground states without magnetic ordering. In this 
paper, we present direct evidence from muon spin relaxation (yz SR) measurements that 
superconducting CeCu 2 .iSt 2 {Tc = 0.7 K) undergoes a random static magnetic ordering 
below T ~ 0.8X. 

It is known that a small amount of off-stoichiometric excess Cm helps to stabilize 

* V 

the superconductivity of CeCu 2 Si 2 . Therefore, we prepared a polycrystal specimen of 
Ce.Cu 2 . 1 Si 2 with the method described in ref. 4. The superconducting transition temper- 
ature Tc = 0.7iC was determined by a resistivity measurement on a small piece cut out 
from the present specimen. A neutron scattering measurement on the crystal structure 
confirmed that there is no minor phase within the accuracy of a few volume percent. 

Zero- and longitudinal-field fj.SH measurements were performed using polarized pos- 
itive muon beams at AGS/BNL and SIN (Zurich) muon channels. In the zero-field mea- 
surements, very small depolarization of muon spins was observed above T = 0.9if, while 
the depolarization rate increased rapidly with decreasing temperature below 0.8if. The 
zero-field muon spin relaxtion functions C 2 {t) (ref. 5) observed at different temperatures 
show no precession signal but have line shapes similar to those observed in dilute-alloy 
spin glasses 5 . This indicates that the magnitude of local fields at muon sites varies widely, 
in contrast, to the cases for uniform ferro- or antiferromagnets. 
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A phenomenological form for the muon spin relaxation function 


<?*(<) = ~exp{~^ch 2 ) + ^exp(-At), ( 1 ) 

with A\/A ~ 2/3 and A^JA ~ 1/3, gives good fits to all the data observed in zero-field. 
The first (second) term of eq. (1) corresponds to the depolarization of muon spins by the 
components of internal fields perpenducular (parallel) to the initial polarization direction 
of the muon spins. Figure 1 shows the temperature dependence of the relaxation rate a. 

In order to distinguish whether this depolarization is due to static or fluctuating 
local fields, we have also performed /zSR measurements by applying longitudinal external 
magnetic fields El = 250 G and lkG parallel to the initial muon spin direction. The muon 
spin polarization had a finite value G z (t) ~ 0.45 with El = 250G, and G z {t) ~ 0.95 
with El — lkG, almost independent on time between 0.5 psec < t < 2 fisec at T = O.liT. 
A similar decoupling of the random field was observed around T = 0.8i£. These results 
indicate that the depolarization shown in Fig. 1 is due predominantly to the static random 
local fields. The rapid increase of a below T ~ 0.8 K then corresponds to the sharp onset of 
static magnetic ordering around the ordering temperature Tm = 0.8iC. Due to the limited 
accuracy of the temperature measurements with the cold-finger dilution cryostat used in 
the present experiment, it is not clear whether the magnetic and superconducting orderings 
occur simultaneously at the same temperature or independently at different temperatures. 
It is, however, evident that the superconductivity and magnetic ordering coexist below 
T c = 0.7 JC. 

The wide distribution of the static random local fields, as observed in the present 
experiment, can be expected either for spin glass (SG) or incommensurate spin-density- 


3 



wave (ISDW) systems. The spin glass ordering is consistent with a large field dependence 
of the susceptibility observed in the non-superconducting system CeCux^St^ (ref.2). If 
one assumes that a majority of Ce (or Cu ) moments participate in the spin-glass freezing, 
the zero-field relaxation rate a ~ 10/isec -1 (the width o/ln ~ 120G of the local field) 
observed at T = 0.05 K corresponds to the dipolar field from a static moment of the order 
of 0.1 Bohr magneton per formula unit. If the small population of the Ce z+ ions with an 
ordered moment of 5 ng forms a spin glass, this value of a is expected for the freezing 
moments on a few percent of the entire Ce atoms. From the present work alone, it is not 
possible to distinguish among the above-mentioned three possible spin structures (the two 
types of SG states and the ISDW state) of CeCu 2 .iSi 2 below Tm- 

Recently, a heavy-fermion superconductor UPtz ( Tc ~ 0.5i£j was found 10 to or- 
der magnetically below T ~ 5 K with an extremely small averaged moment of 0.001 ~ - 
0.02(jtB/U (refs 6,7). URu-iSi-i is another superconducting heavy-fermion system (Tc ~ 
1.0i£) which orders antiferromagnetically below Tm ~ I7if with an ordered moment of 
0.03^a per Uranium atom 8 . With the present results on Ce(7u2.i5t2> we now have three 
heavy-fermion superconductors which show coexisting magnetic ordering with extremely 
small ordered moments (0.001 ~ 0.1 /^b)- This may then be a common feature of the 
superconductivity in heavy- fermion systems. 

This work is supported by the Division of Materials Sciences, US Department of 
Energy under contract 76-AC02-CH00016, the National Science Foundation under DMR 
8503223 and INT 8413978, NASA under NAG-1-416, and by the Japanese Ministry of 
Education, Science and Culture under the Grant-In-Aid for Scientific Research. 
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FIGURE CAPTIONS 


1. Muon spin depolarization rate a, as defined in eq. (l) , derived from the relaxation 
functions observed in CeCu. 2 .iSi 2 in zero field. The onset of magnetic ordering is seen 
around Tm ~ 0.8K. The superconducting transition temperature Tc, determined by 
' a resistivity measurement on a piece cut off from the present specimen, is indicated 
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Motional narrowing of the transverse-field muon spin rota- 
tion signal has been observed in 7 -TiH x for x=1.83, 1.97, and 
1.99. An analysis of the data for TiH^^g near room temperature 
indicates that the mechanism responsible for the motion of the 
muon out of the octahedral site is thermally activated diffusion 
with an attempt frequency comparable to the optical vibrations 
of the lattice. Monte Carlo calculations to simulate the effect 
of muon and proton motion upon the muon field-correlation time 
have been used to interpret the motional narrowing in TiH]_ # g 7 
near 500 K. The interpretation is dependent upon whether Bloem- 
bergen, Purcell, and Pound (BPP) theory or an independent spin- 


pair relaxation model is used to obtain the vacancy jump rate 
from proton NMR T]_ measurements. Use of BPP theory shows that the 
f ield-correlation time can be obtained if the rate of motion of 
the muon with respect to the rate of motion for the protons is 
decreased. An independent spin-pair relaxation model indicates 
that the field-correlation time can be obtained if the rate of 


motion for the nearest-neighbor protons is decreased. 



I 



I. Introduction 


This paper presents an interpretation of the motional nar- 
rowing of the pSR signal in TiH x near room temperature for x = 
1.99 and near 500 K for x = 1.83. We report further analysis of 
the data of Kossler et al. 1 . 

The high mobility of the hydrogen nuclei in metals has at- 
tracted the attention of scientists for over 100 years. 2 During 
the past 35 years, the technique of NMR has been used to deter- 
mine the rate at which these hydrogen nuclei diffuse through 
metals and the mechanism responsible for this motion. 3 To accom- 
pany this, theorists have developed models for diffusion of light 

interstitials 4 . The development of these models has not only been 

) 

helpful in the study of diffusion of hydrogen nuclei, but also 
other light particles such as the muon, which is often considered 
a light isotope of hydrogen (m p « m p /9) due to similarities in 
charge and spin. 

Fukai et al. 5 studied the isotopes of hydrogen in a-NbH x Dy. 
They found that the diffusion coefficient of protons decreased 
and that the proton activation energy approaches that of deuter- 
ons as more protons are replaced by deuterons. The muon extends 
markedly the range of masses studied. Up to now, the muon spin 
rotation (pSR) technique has been used to study ZrH x 6 , VH X 6 , 
NbH x 7 , NbD x 7 , PdH x 8 , TiH x ^ , and YH X 9 . The results of experiments 
on FCC dihydride compounds (TiH x , YH X , and ZrH x ) exhibit the 
following characteristics: partial occupation of interstitial 
(octahedral) sites by the muon at low temperatures, escape from 
these sites and transfer to vacant substitutional (tetrahedral) 



sites in the hydrogen sublattice as the temperature increases, 
and motional narrowing of the /xSR frequency spectrum at higher 
temperatures due to motion of the protons, muons, or both. The 
results of studies on the other hydride compounds indicate that 
the muon. occupies the vacant substitutional sites of the hydrogen 
sublattice upon entering the sample and that motional narrowing 
occurs as the temperature is increased, but it should be men- 
tioned that the vacancy concentration in the hydrogen sublattice 
of these compounds was of the order of 30% or higher. Although 
there exist differences between the /jSR studies of the various 
hydrogen-alloy systems, all of the compounds (except VH 0 . 503 ) 
display an activation energy (as determined from the field-cor- 
relation time at the site of the muon) lower than that measured 
for the proton with NMR (Fig. 1) . This is indeed surprising since 
one naively expects the motion of the muon to be 'impeded by that 
of the proton and that the measured activation would be equal to 
that of the proton. 

In section 2 we present our analysis of the data 1 which 
fields a mean time of stay for a muon at an 0 site. In section 
3 , we develop and show the results of some Monte Carlo calcula- 
tions of the field-correlation time for a muon at a substitution- 
al site as a function of vacancy hopping rate. The results of 
these Monte Carlo simulations are compared to the TiH ^.97 data (T 
» 500 K) of Kossler et al. 1 . 



.II. Activation of Muons Out of Octahedral Sites 

The results of Kossler et al. 1 for pSR studies of TiH x are 
shown in Fig. 2. They find that the muon depolarization rate, a, 
(G x (t)=exp(-A 2 t 2 ) in TiH x for T < 200 K can be parameterized by 


A 2 = Ag(l-P) + PA$ 


( 1 ) 


where a 0 and Ay are the depolarization rates for a muon in oc- 
tahedral and tetrahedral sites, respectively, P is the prob- 
ability of occupying a T site, and 


P = 8(1 - x/2) . 


( 2 ) 


where x is the H/Ti ratio. The term inside the parenthesis is 
equal to the fractional vacancy concentration in the hydrogen 
sublattice and the factor of eight arises from assuming that the 
muon initially stops in an O site and jumps to one of eight 
nearest-neighbor T sites if vacant. The probability of a muon to 
occupy an O site is 1-P, and thus these probabilities are 0.96, 
0.88, 0.32 for x=1.99, 1.97, and 1.83, respectively. In order to 
understand the mechanism responsible for activation of the muon 
out of an octahedral site, let us turn our attention to the 
study of the TiH^gg results since this system has the largest 
probability for a muon to occupy an 0 site. 

The dip in the depolarization rate for TiHj.gg near room 
temperature is similar to that seen in niobium where the depres- 
sion in a was due to muons activating out of relatively more 
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abundant shallow traps, and then finding deeper traps of lower 
concentration. 10 In niobium, the muons in the two traps have the 
same depolarization rate; whereas in TiH^.gg, the muons in the 
abundant and shallow O-site traps have a larger depolarization 
rate than those in the sparser and deep T site traps, yielding an 
asymmetric dip in the depolarization rate. 

The data for those temperature points (330 < T < 370 K) in 
the dip of the depolarization rate (Fig. 2) have been refit using 
a form (Abragamian) for the relaxation function which is better 
suited to motional narrowing than the Gaussian form. 


G x (t) = exp[-A 2 r2{exp(-t/r c ) -1 + t/r c ) ] 


(3) 


where a 2 is the second moment of the field distribution and r c is 
the field-correlation time. Any relaxation function (including 
the one above.) can be represented as: G x (t) = exp[- 7 (t)], which 
allows the definition of an average muon spin decay rate, d 7 /dt. 
When this rate is statistically weighted according to an exponen- 
tial radioactive decay, we can define a dimensionless depolariza- 
tion parameter: 


r 


a = 


exp(— t/r M ) (d 7 /dt) dt, 


(4) 


i J 


where is equal to the mean muon lifetime, 2.2 ^s. For the 
relaxation function listed in Eqn. 3, 


a = A 2 r 2 r c /(r c + r^) . 


( 5 ) 



'a of Eqn. 4 can be naturally related to many theoretical 
descriptions of depolarization. We have found that a is not 
strongly dependent on the form chosen for 7 (t) , and have used 
7 (t) from Eqn. 3 for the purposes of obtaining a through Eqn. 5. 

We used a to obtain the mean time of- stay at a particular site. 

For jiSR studies of Al(Cu), Kossler et al . 11 invoked a 
motional-narrowing theory for 7 (t) to explain the diffusion and 
trapping of muons in terms of correlation functions: 

p* p*' 

7 (t) = Si a|\ dt' dt" Fi(t») Gj_ ( t * — t" ) . ( 6 ) 

'Jo* J 0 

where a|, Fi(t) and Gi(t) are the second moment, probability of 

\ 

occupation and the autocorrelation function for a given site of 
type i, respectively. The summation extends over all interstitial 
site types, which for titanium hydride are the 0' and T sites. To 
define our model we proceed as follows: 

i) We assume that the muon initially stops randomly in 
octahedral sites and then proceeds immediately to a 
near neighbor tetrahedral site if one is vacant, so 
that the initial probability of occupying an octahedral 
site is equal to the probability that it would have no 
near neighbor vacant: (l- 8 c), where c is the vacancy 
concentration . 



ii) Since the results reported by Kossler et al. 1 indicate 
that the muon does not activate out of a T site at room 
temperature, then 

G 2 (t) = 1. (7) 

The autocorrelation function for an 0 site (the shallow 
trap) is 

Gi(t) « exp(-t/r) (8) 

where r is a mean time of stay at an 0 site. 

iii) The transition rate r for a muon to migrate from 0 
sites to a T site is 8c/r, and there is no back migra- 
tion to 0 sites. Thus Fo(t) = 8ce~ rt and F T (t) = (l - 
See” 1 " 11 ) . 

Using these assumptions, we can solve for the mean time of stay 
of the muon at an octahedral site as a function of a by sub- 
stituting Eqn. 6 into Eqn. 4 

r = 8Cr^[-0.5{l + 8c + (l-8c)A^r^/ (a - A^t^) } + 

± [ 0 . 25 { ( 1 - 8 C ) a/ (a - A^2)}2 + 

(1 - 8c)8CAfr2/(a - A^r J) } ] 1/2 ] "1 (9) 

and using the experimental a calculated from Eqn. 5. The ± arises 
because the same a will occur for either slow 0-0 hops and weak 
migration to T sites, or for more rapid hops and greater migra- 
tion. 



We have used the fact that q(T) is single valued at its 
experimental minimum and theoretically when the ± bracket is zero 
to determine the effective vacancy concentration, c. Upon sub- 
stituting in 0.059(3) ps -1 , 0.127(3) /is” 1 , 0.185(3) pS and 2.2 
for a^, aJ, a and r^, we find c to be 0.015(6) which corresponds 
to x=l .97 (1) and agrees with the known concentration of 1.99(2). 

Using the calculated value for the effective vacancy con- 
centration,' we calculated r s (T) and show it in Fig. 3. The line 
is a least squares fit to an Arrhenius form, 

r = r 0 exp(E a /kT) (10) 

where r 0 * io*" 13 ( 1 ) s and E a = 0.48(8) eV. The activation energy 
is a measure of the amount of energy necessary for the muon to 
activate out of an octahedral site and diffuse between o sites. 
The optical vibration frequency i/p for protons in the tetrahedral 
sites of TiH 2 is 3.3 x 10 13 s”*. 12 The muon's oscillation fre- 
quency in the octahedral sites should be close to this value. The 
vibration frequencies for the titanium atoms are about an order 
of magnitude slower. An attempt frequency for hopping should be 
multiplied by a factor representing the number of paths available 
for hopping, a number of order 6. We can therefore only con- 
clude that the mechanism responsible for the activation of the 
muon out of the octahedral site is associated with the vibrations 
of the muon or hydride lattice atoms. 



III. Monte Carlo Simulations 

The disagreement between activation energies measured by pSR 
and NMR has persisted for almost ten years with the one exception 
being the study of VHq.503 * 6 Richter et al . 7 proposed that the 
motion of the muon was limited by the availability of vacancies' 
near the muon in NbH x and NbD x and that the field-correlation 
time, r c , was a measure of the length of time for the field to 
change due to the motion of the muon in the dipolar fields 
generated by the niobium nuclei. The motion of the protons in 
this system primarily has the effect of providing vacancies near 
the site of the muon and will not contribute significantly to the 
field correlation for the muon in this system. However, this is 
not true for the pSR studies of ZrH x 6 , PdH x 3 , TiH x ^ , and YH X 9 
since the dominant contribution to the dipolar field at the site 
of the muon is the sublattice of protons. This means that the 
motion of the protons changes the field at the site of the muon. 
In other words, the motion of the protons alone will produce a 
finite field-correlation time for a stationary muon. Since the 
muon is much lighter, and should be able to jump at a higher rate 
than a proton, the motion of the muon should be impeded by the 
motion of the hydrogen nuclei yielding an activation energy for 
the muon field-correlation time of the same order as that for 
protons. However, since the activation energies and prefactors 
were different for proton- and muon-field-correlation times , 1 we 
implemented Monte Carlo simulations to describe the effects of 
motion of the muon and the protons upon the muon field-correla- 
tion time. 



Since the path for motion of vacancies (hydrogen atoms) in 
TiH x is known to be between nearest-neighbor T sites (as deter- 
mined by NMR), 13 we developed a model with four adjustable para- 
meters: p, p*, Pj_ n , and Pout' to describe the field-correlation 
time for the muon with the restriction that the motion of the 
muon be limited to nearest-neighbor jumping. The motion of the 
vacancies in the unperturbed lattice, i.e. not near a muon, 
occurs with some rate p. When the vacancy is next to the muon, 
the rate for the muon to move or the rate for the vacancy to move 
to the site of the muon is p 1 . This allows the muon jump rate to 
differ from the proton jump rate. Two other rates, p^ n and p Q ut* 
have been added to account for the attraction or repulsion of 
vacancies due to the presence of the muon. These reflect a change 
in vacancy jump rate as the vacancy approaches the muon. p^ n is 
the rate at which a second-nearest-neighbor vacancy to the muon 
will become a nearest neighbor. p ou t is the rate at which a 
nearest-neighbor vacancy to the muon will become a second-nearest 
neighbor. 

The 9x9x9 lattice used for the simulations had periodic 
boundary conditions which allowed particles to diffuse from one 
edge to the opposite edge. Initially, the muon was placed at the 
center of the lattice with seven (= 1%) vacancies randomly dis- 
tributed. The direction of the external magnetic field was picked 
randomly with respect to the coordinate system of the lattice. 

The spins on the remaining lattice locations (excluding the 
positions of the vacancies) were randomly given orientations 
parallel or antiparallel to the field. After this was done, the 



dipolar magnetic field due to the protons in the lattice at the 
site of the muon was calculated. Having done this, the vacancies 
were allowed to move. After each time step, the dipolar magnetic 
field at the site of the muon was computed. The total number of 
time steps was usually 3000-4000, which was large enough to 
determine the field-correlation time at the site of the muon. 

Each simulation produced the dot product of the dipolar magnetic 
field as a function of time with a unit vector in the direction 

A 

of the external field B ex ^ 

B(t) = Sdip(t) • B ext (11) 

The results of each simulation were then multiplied by the 
dot product of the dipolar magnetic field at time t=0 dotted with 
a unit vector in the direction of the external magnetic field 

A A 

(B(0) =Bdip(0) a ®ext) • <B(t)B(0)> , which is used to calculate 
G(t) , was obtained by averaging B(t)B(0) over all simulations 
with the same values for the four parameters: p, p', Pin? and 
Pout • 

G(t) = exp[- 7 2 (t-t» )<B(t' )B(0)>dt» ] ' (12) 

A field-correlation time, r c , is defined such that 
<B(r c ) B(0) >/<B(0) B(0) > = 1/e. 

In order to make the results independent of p, r c was multi- 
plied by p. The results of the simulations are shown in Table I. 
The statistical error in these results is estimated to be of the 
order of five percent or less. 



IV. High Temperature Motional Narrowing of the ^SR Signal 

The results of the Monte Carlo simulations were given in 
terms of a p r c product where p is the vacancy jump rate in an 
unperturbed lattice and r c is the field-correlation time. In 
order to compare the data to the results of the simulations, we 
must multiply the muon field-correlation time by the vacancy jump 
rate. 

The vacancy jump rate can be determined from the proton NMR 
data on 7 -TiH x of Korn and Zamir 12 and of Bustard et al. 13 Korn 
and Zamir 12 , using Bloembergen, Purcell and Pound (BPP) theory, 
found that the proton field-correlation time is: 

V 

T c - (r 0 )/((2-x)/2)exp[E a /kT] s. (13) 

Where r 0 = 2.8(5) * lO" 14 , and E a = 0 . 507 (10) eV. - 
For proton-proton dipolar interactions, the mean dwell time at a 
site, t , (also known as the autocorrelation time) is approximate- 
ly equal to 2 r c . Since 4prc (where c is the vacancy concentration 
which is equal to (2-x)/2 ) is of the order of 1 for an SCC lat- 
tice, then 

p = (8r 0 ) _1 exp[-E a /kT] a" 1 . (14) 

Bustard et al. 13 state that BPP theory may not yield the correct 
values for the field-correlation time, because it is not depen- 
dent upon the type of lattice and does not provide for correla- 
tion effects between diffusing atoms. Instead of using BPP 



theory, they used an independent spin-pair relaxation model to 
interpret their spin-lattice relaxation times, T]_. Their Monte 
Carlo model monitored the time evolution of a large number of 
spin pairs. Using the results of their simulations, they deter- 
mined the mean dwell time for a proton at a site, r, from their 
data and also from the data of Korn and Zamir 12 and found 

r = (9.3(6) * 10“ 15 ) / ( 1/2 (2-x) ) exp [0.555(7) eV/kT] s (15) 

This yields a vacancy jump rate of 

p= (2.7(2) * 10 13 ) exp[-0. 555 (7) eV/kT] s" 1 . (16) 

Multiplying the vacancy jump rates obtained from both references 
by the muon field-correlation times for TiHi. 97 , we obtain the 
p r c product values shown in Table II. The proton hopping rates 
obtained by the later and more sophisticated work of Bustard et 
al . 13 are roughly half those of Korn and Zamir 14 . 

The pt- c products obtained using either set of proton hopping 
rates are only constant below 500 K and rise to about twice their 
lower temperature values by 563 K. Naively, one would expect the 
vacancy hopping rates near the muon to be about equal to the rate 
further away: p ou t “Pin = P ; and the rouon to hop to a nearest 
neighbor vacancy at least as fast as a proton: p' > p. We see 
from a comparison of the tables this does not work for T > 500K 
for either set of proton hopping rates^and for those of Bustard 
et al. (Table II, column 2) at any temperature. An increase of 



p r c , but with p 1 > p can be achieved by a reduction of the rate 
at which a vacancy hops to a nearest neighbor site, simulation 8 
or 10, pj_ n = p/10, otherwise one must also reduce p', simulation 
7 . A very dramatic increase in the rate of vacancy hopping away 
from a nearest neighbor site can increase p r c moderately, simula- 
tion 5. However, if the results of Bustard et al. are used, we 
see that reducing the muon's motion alone (p' = 0, simulation 7) 
is not sufficient and the approach of a vacancy must be 
restricted, pj_ n < p, see simulations 7-12. The treatment by 
Bustard et al. was developed specifically for motion of the 
hydrogen nuclear spins on an SC lattice and would appear to be 
better suited for extraction of a proton autocorrelation time 

* 

from T]_ data compared to BPP theory. However, using an indepen- 
dent spin-pair relaxation model to obtain the vacancy jump rate 
requires that the muon affect the motion of the nearest-neighbor 
protons . 

A vacancy to muon-site hopping rate less than that for a 
proton implies the involvement of some other mechanism. This may 
also be related to the requirement that p-j_ n be reduced relative 
to p. If the nearest-neighbor hydrogens are bound to the muon by 
an excess of 0.1 eV this would reduce Pi n sufficiently and pos- 
sibly force muon hopping to next-nearest-neighbor sites, which 
are further apart and would hence reduce the prefactor by an 
order of magnitude. 



V. 


Conclusions 


We have reanalyzed the TiH]_ # gg data taken near room tempera- 
ture by Kossler et al . 1 and have shown that they can be fit to a 
model with many shallow traps and a low concentration of deeper 
traps for the muon. We find that the mechanism responsible for 
the activation of the muon out of the shallow trap (O site) is 
associated with the vibration of the -hydrogen sublattice. Second- 
ly, a method of simulating the effect of the motion of protons 
and the muon on the field-correlation time is presented and 
applied to the motional narrowing of the ^SR signal in TiH 1>97 
for temperatures near 500 K. The conclusions from the comparison 
of the p r c obtained from the data and from the simulations are 
dependent on the method of analysis to obtain an autocorrelation 
time for the proton from proton NMR data for 7-TiH x . Using the 
results of Korn and Zamir 12 , we find that the motion of the muon 
decreases with respect to the motion of the proton as the temper- 
ature increases. Using the results of Bustard et al. 13 , we find 
that a stationary muon does not explain the results of Kossler et 
al . 1 and that the presence of the muon inhibits the motion of the 
nearest-neighbor protons. It would seem that the method used by 
Bustard et al . 13 is better suited for analysis of proton data 
for 7-TiH x , but the model needed to describe the results of nSR 
experiments is more complicated than the model developed using 
BPP theory employed by Korn and Zamir 12 to derive the vacancy 
jump rate. We believe that the conclusions reached for TiH^.gy 
are also valid for TiHi # 33, since the proton and muon correlation 
times scale inversely with the vacancy concentration. We -also 



believe that these conclusions are valid for other FCC metal 
hydrides, i.e., ZrH x and YH X . 
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TABLE I prg products from Monte Carlo simulations, p is the 
vacancy hopping rate far from the muon, p* is the rate at which 
the muon moves to a nearest-neighbor vacancy, pj_ n is the vacancy 
rate for next-nearest-neighbor to nearest-neighbor jumps, p ou t is 
the rate for the converse process of Pi n , and r c is the field 
correlation time for the muon. 


Simulation# Probabilities or c 

01 P'=Pout=P' Pin =1 °P 6 

02 Pin~Pout = P' P ,= 1°P 12 

03 P’=Pin=P/ Pout=P/100 13 

04 P ,= Pin~Pout = P 23 

05 P ,= Pin~P» Pout =1 °P 32 

06 Pin"Pout = P' P* “P/100 49 

07 Pin“Pout~P» P ,= 0 52 

08 Pout=P/ Pin = P/l°/ P'“10p 61 

09 Pin"P' P out = P/ 100, p'=0 61 

1° Pout=P’ = P' Pin“P/l° 98 

11 Pout=P'=P/ Pin=P/100 450 

12 Pout=P/ Pin = P/l°°/ P'=0 452 



TABLE II. p r c products for TiH 1<97 data. 


T (K) 

Tq (uS) 

4.3(7) 

P (K&Z ) T ri 

P(B) Tr> 

432 

24(9) 

39(10) 

458 

2.1(3) 

25(9) 

44(11) 

477 

1.3(2) 

26(9) 

49(11) 

493 

0.9(1) 

26(9) 

50(11) 

508 

0.8(1) 

34(11) 

68(15) 

523 

0.8(1) 

45(14) 

93 (20) 

543 

0.49(9) 

43(15) 

93(22) 

563 

0.35(7) 

45(16) 

101(26) 

583 

0.04(6) 

7(H) 

16(26) 


K&Z indicates that the vacancy hopping rate has been taken from 
Korn and Zamir 12 . 

B indicates that the rate is from Bustard et al. 13 . 
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Abstract 

A search for nuclear-bound states of the n meson has been carried out. Targets 
of lithium, carbon, oxygen, and aluminum were placed in a a + beam at 800 MeV/c. 
A predicted . rj bound state in n ( E x ~ 540 MeV) with a width of «9 Me".' was 
not observed. A bound state of a size 1/3 of the predicted cross section would 
have been seen in this experiment at a confidence level of 3a ( P>0 .998 7 ) . 


* Now at Osaka University, .Japan 
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This paper describes a search for a novel nuclear excitation involving the 
creation of a bound n-meson in the nuclear medium. The concept is similar in 
spirit to a number of ideas which have been recently vigorously pursued. Some 
familiar examples are A hypermiclear states, F. -hypernucl ear states, antiprotonic 
nuclear states, and various diharyon resonances. In each case an attractive 
particle-nucleus potential is required and some mechanism to inhibit the decay 
process, such as strangeness conservation in the case of the A. 

Several suggestions of the existence of bound states of the q-meson in a 
wide range of nuclei have recently been published 1-3 . The suggestions of this 
novel nuclear excitation are based on bound state formation through the attrac- 
five N-q channel of the N (1535), where N' (.1535) is the (uN) resonance with 
(I.J 11 ) = (l/2,l/2“) and mass 1535 MeV/ c . This resonance dominates q pro'duc- 
tion near threshold. Bhalerao and Liu have shown, by a coupled-channel analy- 
sis, that the low-energy qN interaction is attractive with a scattering length 
of (0.28+0. 202i) fermi. The attractive interaction is a consequence of the 

x 

threshold being below the '’N(1535) resonance. 

Liu and his collaborators have examined the consequences of this attractive 
interaction in the formation of a bound-p state as a function of mass number. 
Their study indicates that nuclear bound states could exist for mass numbers 
larger than A=10. At low mass numbers, only s-state bound p’s are predicted. 
At larger mass numbers, p- and d-states could become bound. Both binding 

energies and widths increase with A. The optimum case, in their analysis, is 
l "0, formed from the (x + ,p) reaction on 1(3 0 at a momentum near 740 MeV/c." 1 

At an angle near 15°, the momentum transfer is favorable for the transition 
involving the conversion of p-shell neutron to an s-shell p. For higher mass 
numbers, the increase in predicted width would make this excitation more 
difficult to see over the continuum (c + ,p) background which is present. 



An experiment to test these predictions was devised with the positive pion 
heam available at the Low-F,nergy Separated Beam I at the Brookhaven \1.-te rnat Lng 
Gradient Synchrotron (AGS), and the _ Moby Dic.k Spectrometer. The experimental 
arrangement is virtually identical with that used for the production and 
measurement of hyperuuclei, and it has been described in detail in a number of 
publications (see Ref. 5 and related references). The only differences involve 
the selection of pions , rather than kaons , in Che incident particle channel and 
protons in the exiting particle channel. 


The 

spectrometer was 

set 

at 15° to 

be 

near the 

maximum for bound 

eta pro- 

duction 

as predicted in 

the 

analysis 

of 

Liu and 

1 1 

Haider >' . For 

example, 


reference 3 predicts the production of an p excitation of width of 9 MeV ( F*WHM) 
and a peak cross section of about 30 ub./sr-MeV, for an incident pion momentum of 
300 MeV/c on an 16 0 target. The peak would appear near zero binding energy in 
the (^O+p) missing mass spectrum corresponding to the emission of protons with 
248 MeV kinetic energy in the lab frame, and an excitation of 540 MeV in the t5 0 
system. 

Four targets were selected for examination; the target parameters are 
listed in Table I. The oxygen target was in the form of water. For lithium, no 
p-bound states are predicted, while for carbon the binding is predicted to be 
marginal. Oxygen is expected to display the. largest bound state cross section; 
for larger A the cross sections are smaller, while the widths grow larger. 

i 

To confirm a cross section scale and to measure a spectrometer momentum 
acceptance function, the reaction p(r’ ! ’,p)it + was measured for a :c' r momentum of 
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525 MeV/c. This reaction has been measured by the Leningrad group of Gordeev et 
al. 8 , and their reported cross sections were in reasonable agreement (-20%) with 
our measurements. To establish the acceptance function, the nominal momentum 
setting of the proton spectrometer was varied from 520 to 780 MeV/c, and the 
strength of the observed n + missing mass peak was used to determine the 
acceptance. It was desirable to limit the spectrometer acceptance correction to 
no more than 30% of the central value over the entire spectrum. With this 
criterion, an acceptance range of =80 MeV/c was obtainable; i.e., the relative 
acceptance is .everywhere larger than 0.7 throughout the range. To obtain a 
sufficiently broad range in outgoing proton momenta, overlapping runs were taken 
for spectrometer central values of 657, 700, and 740 MeV/c. 

The spectrometer resolution was measured by analyzing (p,p') events' for 


12 


C, recorded simultaneously with the (r + ,p) events. The missing mass resolu- 


tion was studied in two separate ways: by using TRANSPORT matrix elements to 
calculate the particle momenta, and also by using the program RAYTRACE 8 , which 
includes, magnetic corrections to all orders. The spectrometer resolution was 
measured to be 4 MeV ( FWHM) using the TRANSPORT analysis and 2.5 MeV (FWHM) 
using RAYTRACE. For either mode of analysis, the resolution is sufficient to 
detect the predicted n-states without serious resolution broadening of the peak. 

These (n,p') studies also serve to confirm the energy scale and the energy 
losses in the target and beam windows of the experiment. 

The results of the experiment on the 4 targets are shown in the composite 
diagram of Fig. 1. The inclusive proton spectra of this experiment show a 
qualitative similarity to recent (p,p) reaction studies on various nuclei by 
Garreta et al. We know of no comparable (r ,p) data at these pion energies. 
The data for each target show a smooth (tu + ,p) cross section down to an energy 
corresponding to the q production threshold. At lower energies, the cross 
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525 MeV/c. This reaction has been measured by the Leningrad group of Gordeev et 
al. 6 , and their reported cross sections were in reasonable agreement (=20%) with 
our measurements. To establish the acceptance function, the nominal momentum 
setting of the proton spectrometer was varied from fS 20 to 780 MeV/c, and the 
strength of the observed rc + missing mass peak was used to determine the 
acceptance. It was desirable to Limit the spectrometer acceptance correction to 
no more than 30% of the central value over the entire spectrum. With this 
criterion, an acceptance range of =80 MeV/c was obtainable; i.e., the relative 
acceptance is everywhere larger than 0.7 throughout the range. To obtain a 
sufficiently broad range in outgoing proton momenta, overlapping runs were taken 
for spectrometer central values of 657, 700, and 740 MeV/c. 

The spectrometer resolution was measured by analyzing (p,p') events' for 
1 2 C, recorded simultaneously with the (Tt + ,p) events. The missing mass resolu- 
tion was studied in two separate ways: by using TRANSPORT 7 matrix elements to 
calculate the particle momenta, and also by using the program RAYTRACE 8 , which 
includes magnetic corrections to all orders. The spectrometer resolution was 
measured to be 4 MeV ( FWH.M) using the TRANSPORT analysis and 2.5 MeV ( FWHM) 
using RAYTRACE. For either mode of analysis, the resolution is sufficient to 
detect the predicted n-states without serious resolution broadening of the peak. 

These (p,p’) studies also serve to confirm the energy scale and the energv 
losses in the target and beam windows of the exoeriment. 

The results of the experiment on the 4 targets are shown in the composite 
diagram of Fig. 1. The inclusive proton sneer ra of this experiment show a 
qualitative similarity to recent (p,p) reaction studies on various nuclei by 
Garreta et ai. We know of no comparable (~ ,p) data at these pion energies. 
The data for each target show a smooth (m + ,p) cross section down to an energy 


corresponding to the q production threshold. 


At Lower energies. 


the cross 
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section appears to increase, and it is plausihl.e to attribute th is 'increase to 
the quasi-free q production process. For the aluminum, oxygen, and possibly 
carbon cases an eta peak would be expected to appear near the position of the 
arrow. 

The sensitivity of this experiment to narrow q bound-state peaks is 
obviously compromised by the large proton continuum background observed in this 
experiment. This background is presumably attributable to nuclear protons 
ejected by the incident pions through various processes, including quasi-Cree 
knockout, multiple pion and proton scattering, and pion absorption. To estab- 
lish the experimental sensitivity quantitatively, a statistical analysis was 
carried out in detail. The analysis was made on data that were not corrected 
for momentum acceptance, since the correction process increased the spread of 
the data points significantly. Fluctuations about a polynomial fit to the 
uncorrected data were analyzed with a standard least squares fitting code. 

The quantitative statement of the experimental sensitivity, based on the 
largest observed fluctuation in the data,, is the following: the detection of a 
peak with a full-width at half maximum of 9 MeV in the ^0 data occurs at a 
confidence level of 3 cr (0,9987) for a peak height of 8.7 ub/sr/MeV. This size 
is about one-third of the prediction of Ref. 3. The ^0 target was predicted in 
Ref. 3 to be the most favorable case, i.e., the one displaying the largest cross 
section. Similar sensitivities obtain for the lithium and carbon cases, while 
for aluminum the size of the fluctuations, due to a poorer statistical accuracy, 
preclude any strong statements. 

It is interesting to remark on the similarities of the spectra shown in 

Q 

Fig. 1 with the spectra obtained in the experiment of Garreta — both of which 
experiments were designed to search for narrow states located near the onset of 
a quasi-free process. In the latter the search is for p states, while in our 
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work, we are searching for a structure due to hound q states. 

The (p,p) process is thought to proceed via the annihilation of the p with 
a target nucleon to produce, on average, 5 pions which subsequently interact 
with the A- 1 target nucleons and eject protons. Thus the subsequent stages of 
the process are similar to the (rc + ,p) reaction. Reference 9 documents the fact 
that over a very wide range of proton momenta, the cross section has a 
Maxwellian shape: 


d 2 a/dQdE = C( E) 1/2 exp{ -E/T} 


where E represents proton kinetic energy. 

We adopted the same parameterization to characterize the (iT + ,p) data. ' For 
the fitting procedure, it was necessary to exclude the q quasi-free region; 
hence only data outside the allowed kinematical range of q production was used. 
The resulting fits were then extrapolated into the q quasi-free region, these 
fits are indicated in Fig. 1. We believe that the obvious excess which occurs 
near the onset of the q threshold is attributed to q quasi-free processes. 

An expanded region near the q threshold is shown in Fig. 2 for the oxygen 
target. * To indicate the various reaction mechanisms more explicitly, three 
curves are shown: the dashed curve is the Maxwellian fit to the (x + ,p) inclu- 
sive process, the solid curve shows an estimate of the nuasi-free q production 
added to the Maxwellian, and the docted curve shows the p-shell bound eta state 
predicted by Liu and Haider »' . It should be clear from the figure that a 
peak of the predicted size would be visible in the experiment; it is not seen. 

The size of the quasi-free production observed in the ^0 data is roughly 
in Line with the observations of Peng and collaborators'^ at a somewhat lower 
heam momentum of 680 MeV/c. An estimate of the magnitude of the expected 



ORIGINAL PAGE IS 
OJE EOOR QUALITY. 

quasi— f ree eta production was based on 1 1 10 elementary cross sections of 
l l l o 

Litchfield et al. and of Brown et al. “ (the latter for the charge conjugate 
(n - ,n) reaction). The shape of the onset of the quasi-froe process near 

threshold was taken from Ref. 1, and that result was arbitrarily normalised to 
our data as shown in Fig. 1. From this normalization, a result for the 

effective number of scatterers, defined as, 

,, do(nuclear) 

eff = 

do (elementary) 

can be obtained. The preliminary estimate of the quasi-free cross section, 

derived from the fit of Fig. 2 and integrated over the spectrum shape presented 
by Ref. 3 leads to a value do/dQ (qf) = 150 ph as compared to a value 191 pb 
derived from Ref. 11. We estimate therefore that N e ff is of the order of 
unity (=0.8). This is quite comparable to the value characteristic of hypernu- 
clear cross sections and an order of magnitude higher than N e ff deduced from 
(p,p) nuclear interactions. 

In summary, the search for a narrow q -nuclear bound state has produced 

negative results. A peak one-third the size predicted would have been seen in 
oxygen at a confidence level of 3o . The conclusion is either that widths for 
such excitations are much larger than predicted, or that the production cross 
sections are much smaller than predicted, or both. 
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1. The proton kinetic energy spectra obtained for the targets exam i.ned in this 
experiment at 800 MeV/c i.nc ident it' 1 '. Tn ear.li case , a Max we 1 Man function 
was fit to the (tt + ,p) inclusive proton energy above the eta production 
threshold. The arrows indicate that threshold for each target. The error 
bar shown indicates the standard deviation for a typical datum near 
threshold . 

2. The region of the oxygen spectrum in which an n-bound state would appear is 
shown in this figure. The dashed curve shows the extrapolated Maxwellian 
fit to the inclusive (>c,p) background, and the solid curve shows the 
quasi-free q-production added to the (t,d) background. The dotted durve 
indicates a state in which a p -^/2 shell, neutron has been transformed into 
a bound rj , with the si a e , width and binding predicted by ref. 2. The data 
are clearly inconsistent with that prediction. 
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Table l. 

I r rad iat ions 


The target thicknesses and the total dead-time-corrected pion 
for Che spent ra obtained in this experiment are indicated. 


Target 


Th ickness 


(gin /cm ' ) 


Dead -Time-Corrected 
Pion Irradiation 


Li 

2.25 

C 

l . 5° 

Water 

2.15 

AJ7. 

2.03 

Polyethylene 

2.15 


( calibration) 


2.92x10 
2.99x10 
A. 57x10 
1.66x10 
1.40x10 


10 

10 

10 

10 

10 
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